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Abstract 



The mathematical methods that have been used to analyze the 
statistical properties of boson fields, and in particular the coher- 
ence of photons in quantum optics, have their counterparts for 
Fermi fields. The coherent states, the displacement operators, 
the P-representation, and the other operator expansions all pos- 
sess surprisingly close fermionic analogues. These methods for 
describing the statistical properties of fermions are based upon 
a practical calculus of anti-commuting variables. They are used 
to calculate correlation functions and counting distributions for 
general systems of fermions. 
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1 Introduction 



The Pauli exclusion principle plays an essential role in describing the behav- 
ior of the particles, both simple and complex, that we now call fermions. It 
is known to play a key role in determining the structure of the most fun- 
damental elements of matter. These are systems like atoms, in which the 
phase-space density of fermions, electrons in this case, is quite high. But 
when fermionic atoms move freely in space or even when they are trapped 
electromagnetically, their phase-space density is usually so low that the ef- 
fects of the exclusion principle remain completely hidden. A number of re- 
cent developments, however, point to the possibility of achieving much higher 
densities of fermionic atoms both in electromagnetic traps and in free space. 

The various methods of optical cooling that have been developed for 
atomic beams work as well for fermions as they do for bosons and produce 
beams with temperatures of the order of 100 fi°K. Cooling fermions evapora- 
tively to still lower temperatures poses a problem that requires a less direct 
solution. Evaporative cooling becomes inefficient for fermions since the ex- 
clusion principle tends to suppress collisions of identical atoms. It may be 
implemented nonetheless by sympathetic means Jl|], e.g., by cooling bosonic 
atoms at the same time, so that energy exchange still takes place freely. 
It seems possible thus that the realization of degenerate Fermi gases may 
become an important byproduct of Bose-Einstein condensation. 

The detection methods that will be used in measurements on beams of 
cold fermionic atoms will be essentially the same as those now used on bosonic 
atoms cooled by optical or evaporative means. The measurements on bosons 
can be most conveniently described, in fact, by mathematical methods that 
were introduced in the context of quantum optics [TJ, |3J. 

Much of the work in quantum optics, we may recall, is couched in the 
language of coherent states, which are eigenstates of the photon annihilation 
operators. They contain an intrinsically indefinite number of quanta but can 
nonetheless be used as a basis for describing all states of the electromagnetic 
field. While pure coherent states are not physically attainable in bosonic 
systems with fixed numbers of particles, it likewise remains useful to describe 
boson fields in terms of suitably weighted superpositions and mixtures of 
coherent states. The weight functions associated with these combinations 
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may be regarded as quasi-probability densities in the spaces of coherent- 
state amplitudes. The function P in the coherent-state representation of 
the density operator |2|, [| plays this role; other quasi-probability densities 
including the Wigner function || [|] and the Q function || || play similarly 
convenient roles in representing the density operator. 

In the case of fermion fields, the vacuum state is the only physically 
realizable eigenstate of the annihilation operators. It is possible, however, to 
define such eigenstates in a formal way and to put them to many of the same 
analytical uses as are made of the bosonic coherent states. Since fermion field 
variables anti-commute, their eigenvalues must, as noted by Schwinger H, 
be anti-commuting numbers. Such numbers are Grassmann variables. They 
can be handled by means of the simple rules of Grassmann algebra J7], which 
we include here so that the calculations may be self-contained. 

Within this context we formulate ways of expressing and evaluating a 
broad range of the correlation functions that are measured in experiments 
involving the counting of fermions. Central to this task is the expression of 
the quantum-mechanical density operator in terms of Grassmann variables. 
We develop a number of ways of doing that in general terms and present 
a detailed discussion of the density operators for chaotically excited fields. 
Included among the latter is a particularly useful Gaussian representation 
of the grand-canonical density operator for fermion fields. Having evaluated 
the statistically averaged correlation functions, we apply them to fermion 
counting experiments and illustrate their use in determining the counting 
distributions. 

We find throughout this work that notwithstanding great mathematical 
differences, many close parallels can be established between the expressions 
evaluated for fermion fields and the more familiar ones for boson fields. In 
particular, for example, we can construct a family of quasi-probability den- 
sities, as functions of the Grassmann variables, with properties parallel to 
those of the entire family of quasi-probability densities for bosons, including 
the P, Q, and Wigner functions. We can then evaluate the mean values of 
ordered products of fermion creation and annihilation operators by perform- 
ing integrations over the Grassmann variables while using the appropriate 
quasi-probability density as a weight function. In both cases, we trade an in- 
homogeneous commutation relation and an ordering rule for a homogeneous 
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commutation relation and a quasi-probability density. For boson fields the 
integrations are taken over commuting variables, which may be treated as 
if they were classical variables. For fermions, on the other hand, the inte- 
grations are over anti-commuting variables that have no classical analogues. 
The weight functions for these integrations are nevertheless in one-to-one 
correspondence with the quasi-probability densities for bosons, so it seems 
appropriate to give them similar names. We have followed that convention 
for several other parallels as well. 



2 Notation 



Let us consider a system of fermions which may be described by the cre- 
ation a* and annihilation a m operators which satisfy the familiar but ever 
mysterious relations 



{(Z n , O m } $nm (1) 

{a n ,a m } = (2) 

{4,4J = (3) 

Gn|0> = 0. (4) 



in which |0) is the vacuum state. 

We shall use lower-case Greek letters to denote Grassmann variables. 
These anti-commuting numbers j n and their complex conjugates 7* satisfy 
the convenient relations 

{7n,7m} = (5) 

K,7 m } = (6) 
{7:, 7; t } = 0. (7) 

We shall also assume that Grassmann variables anti-commute with fermionic 
operators 

{7n, a m} = (8) 

and commute with bosonic operators. And we make the arbitrary choice that 
hermitian conjugation reverses the order of all fermionic quantities, both the 
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operators and the Grassmann numbers. Thus for instance we have 

(a^al^y = -f 4 a 3 p;a\. (9) 

3 Coherent States for Fermions 
3.1 Displacement Operators 

For any set 7 = {7^} of Grassmann variables, let us define the unitary 
displacement operator D(j) as the exponential 

D{l) = exp (ah ~ 7>i) j • (10) 

One of the useful properties of Grassmann numbers is that when, as in the 
preceding definition, they multiply fermionic annihilation or creation op- 
erators, their anti-commutativity cancels that of the operators. Thus the 
operators a J 7$ and 7*Oj simply commute for i ^ j. So we may rewrite the 
displacement operator as the product 

Din) = n ex p(°H-^) ( n ) 

i 

= n[ i+a ^-^+( a ^-i)^]- ( i2 ) 

i 

By the same token, the annihilation operator a n commutes with all the op- 
erators a'ni and 7*0, when n 7^ i, and so we may compute the displaced 
annihilation operator by ignoring all modes but the nth: 

D\j)a n D(~/) = JJexp (7*0* - a n JJexp (a^j ~ Ipj) 

i 3 

= exp (7*a„ - a\n n ) a n exp (a*7 n - 7*a n ) 

= (1 - a* 7n - §7X47n) a n (l + a\n n ~ |4 7n7>n) 

= (l - «n7n - \lnln) d n (l + a+7 n + ±7n7n) 

= a n - a\n n a n + a„a] l 7 n = a n + 7„. (13) 
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Similarly 

Dt( 7 )ai£>( 7 ) = a t+'yj;. (14) 
We may use the Baker-Hausdorff identity 

e A + B =e A e B e -l[A,B]^ (15) 

which holds whenever the commutator [A, B] commutes with both A and B, 
to write the displacement operator D(ct) in forms that are normally ordered 

exp ^XX 7i j exp (^~^2^i a ^j = ex P (^2 ( a hi-l* a i)^j e^*' 7 

Av(7) = £(7)exp(±£ i7 * 7i ) (16) 

and anti-normally ordered 



ex P I _ 5Z 7 '* a y exp ^XM 7i J = exp [Yl ( a * 7i ~ 7 ** ai ) I e *' r ' 7 

D A ( 7 ) = Din) exp(-f£. T * 7 ) , (17) 
in which we have employed the concise notation 

i 

an abbreviation which we shall use occasionally but not exclusively. The 
identity fll5 ) also allows one to show that the displacement operators form a 



ray representation of the additive group of Grassmann numbers, 

D(a) D((3) = D(<x + (3) exp [§ £ 4 - • (19) 

3.2 Coherent States 

For any set 7 = {72} of Grassmann numbers, we define the normalized co- 
herent state 1 7) as the displaced vacuum state 

| 7 ) = £>( 7 )|0). (20) 
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By using the displacement relation (|T3"D, we may show that the coherent state 
is an eigenstate of every annihilation operator a n : 

a n | 7 > = a n D(7)|0) = D(7)Dt( 7 )a n D(7)|0) 

= D(j) (a n + 7n )|0)= J D( 7 )7n|0)=7n J D(7)|0) 

= 7n|7>- (21) 

By using the product formula flT2] ) for the displacement operator, we may 
write the coherent state in the form 



|7> = J D(7)|0) = n[ 1 + «l 

i 

= n( 1 + a ^-lT*Ti)iO) 



7 - -,/"; 1 ( 4 a i ~ \ ) lilt 



|0> 



exp 



2 7i 7? 



|0>. 



(22) 



It may be worth emphasizing that in this formula the creation operator a] 
stands to the left of the Grassmann number 7$. Apart from these ordering 
considerations, this formula takes a form closely analogous to the one that 
defines bosonic coherent states. 



The adjoint of the coherent state I7) is 



( 7 | = <0|£>t( 7 ) = <0| exp [J2 W°< 



|7*7i) 



and it obeys the relation 

(7l4 = (7l7n- 

The inner product of two coherent states is 

( 7 |/3) = exp ( J2 (7*A - |(7*7, + ml) 



(23) 
(24) 

(25) 



so that 



(0l7><7l/3> = exp 



E(A*-7*)(A-7.) 



(26) 



In contrast to the case of bosons, we may for fermions define for any 
set a. = {cKj} of Grassmann numbers the normalized eigenstate \a)' of the 
fermion creation operators a\ as the displaced state 

\<*)' = D(a)\l) (27) 

where |1) is the state in which every mode is filled: 

\l)=l[ai\0). (28) 

n 

By using the displacement relation (|H), we may show that the state \et)' is 
an eigenstate of every creation operator a*: 

aj»' = alD( a )\l) = D( a ) D^( a )aiD( a )\l) 

= D(a) (al + a* n )\l)=D(a)a* n \l) = a* n D(a)\l) 
= <!«)'• (29) 
The adjoint relation is 

\ct\a n = \cx\a n . (30) 
An explicit formula for the eigenstate |a)' follows from its definition (p7|): 

|a)' = Y[(l-a*ai + ±a*ai)\l}. (31) 



3.3 Intrinsic Descriptions of Fermionic States 

The occupation-number description of states of fermions has well-known am- 
biguities. For n m, for example, the state |l n l m ) may be interpreted as 

4 a J»l°) or as a m4lo) = -44jo). 

The creation operators themselves provide an unambiguous description 
of fermionic states, 

\^) = J2c(n 1 ,n 2 ,...)ai i ai 2 ...al m \0), (32) 
W 

which transfers to the coherent-state representation 

(a|<0) = exp (-§ J^n «X) Yl °( ni ' ™ 2 ' • • a «i a "2 • • • a n m ( 33 ) 

{n} 
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without any ambiguity or extra minus signs. Because coherent states are 
defined in terms of bilinear forms in anti-commuting variables, there is no 
need to adopt a standard ordering of the modes. 



4 Grassmann Calculus 



4.1 Differentiation 

Since the square of any Grassmann variable vanishes, the most general func- 
tion /(£) of a single anti-commuting variable £ is linear in £ 

f{Z) = u + &. (34) 

We define the left derivative of the function /(£) with respect to the Grass- 
mann variable £ as 

— = t. (35) 

Note that if the variable t is anti-commuting, then we may also write the 
function /(£) in the form 

f(Z) = u-t£. (36) 
Now to form the left derivative, we first move £ past t, picking up a minus 



sign and obtaining the form fl3~4 ) and the result (^). In this case, the right 
derivative is —t. In the present work, we shall use left derivatives exclusively 
and shall refer to them simply as derivatives. 



4.2 Even and Odd Functions 

It is useful to distinguish between functions that commute with Grassmann 
variables and ones that do not. We shall say that a function f(ct) that 
commutes with Grassmann variables is even and that a function f(ct) that 
anti-commutes with Grassmann variables is odd. We shall often note the 
evenness or oddness of the functions we introduce. 
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4.3 Product Rule 

To compute the derivative of the product of two functions /(<*) and g(oc) 
with respect to a particular variable «j, one may explicitly move the «j in 
g(cx) through f{ot) or one may move the operator representing differentiation 
through the function f(ot). In either case if the function f(ot) is odd, then 
one picks up a minus sign. The product rule is thus 

{f(a)g(a)) = — g(a) + cr(f) f(a) — , (37) 



doti dati doti 

where the sign a(f) of f(oc) is —1 if f(ct) is an odd function and +1 if f(ct) 
is even. 



4.4 Integration 

We define a sort of integration over the complex Grassmann variables by the 
following rules 



da n = J da* n = (38) 
da n a m = 5 nm (39) 
da* n a* m = 5 nm . (40) 

This integration due to Berezin is exactly equivalent to left differentiation. 

We shall typically be concerned with pairs of anti-commuting variables 
«i and a*, and for such pairs we shall adhere to the notation 



d a n = J da* n da n (41) 

in which the differential of the conjugated variable comes first. Note that 

da n da* n = —da* n da n . (42) 
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We have been using boldface type to denote sets of Grassmann variables; we 
shall extend that use to write multiple integrals over such sets in the succinct 
form 

Jd 2 a = J Y[d 2 ai . (43) 

We shall also occasionally employ the concise notation 

a*-/3 = J>;/3 n (44) 

n 

for sums of simple products over all the modes of the system. 
The simple integral formula 

J d 2 a n e~ a p aq = J da* n da n (l — a*a q ) = — J da* n da n a* p a q 

= J da* n a* da n a q = J da* n a* p 5 nq = S np 5 nq (45) 

provides a useful example of Grassmann integration. We also note the general 
rule 

jd 2 af{\a) = \\\ 2 jd 2 (3f{(3), (46) 

in which A is an arbitrary complex number and in which f(Xa) is an abbrevia- 
tion for a function which necessarily depends on both \a and X*a*. This rule 
owes its strange appearance to the definition of integration as differentiation. 

Some further examples are the integral of the exponential function 

J d 2 a exp (P*a + a*7 + aa*) = exp (47) 

and the Fourier transform of a gaussian 

J d\ exp « - + XiC) = A exp (48) 

where A is an arbitrary complex number. The latter integral can be written 
in a somewhat more-general form which is no longer a Fourier transform: 



/ 
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4.5 Integration by Parts 



Let us first observe that the integral of a derivative vanishes 

rf 2 «^ = (50) 
don 

because the derivative with respect to the variable q.{ lacks the variable a*. 
In particular the integral of the derivative of the product of two functions 
also vanishes, and so by using the product rule (0), we have 



d 

d 2 ct — (f(ac)g(cx)) = \ d' 2 ct 



d JM) g(a)+ a U)!(a) a -P 

OOLj I OOCi 



which is the formula for integration by parts 
'df(a 



da 



dou 



g(a) = -a(f) Jd 2 af(a) 



dan 



= 0, 
(51) 

(52) 



where the sign a(f) is +1 if the function /(«) is even and —1 if it is odd. 



4.6 Completeness of the Coherent States 

We may use our Grassmann calculus to show that the coherent states are 
complete. Let us consider the state 

|/) = (c+da+)|0>, (53) 

which for arbitrary complex numbers c and d is an arbitrary single-mode 
state. Then its inner product (7I/) with the coherent state I7) is the correct 
weight function for the coherent-state expansion since 

Jd 2 l(l\f)\l)= /d 2 T ( C + d 7 *) (l + 77*-7« t )|0) = (c + dat)|0) = |/). 

(54) 

The reader may generalize this example to the multi-mode case. The coherent 
states in fact are over-complete. 
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4.7 Completeness of the Displacement Operators 



For a single mode, the identity operator / and the traceless operators a, , 
and \ — a) a form a complete set of operators. Since by using the expression 
( |T2] ) and our Grassmann calculus, we may write each of these operators as 
an integral over the displacement operators 



/ = 


J d 2 7 77*D( 7 ) 


(55) 


a = 


y"d 2 7 (-7)^(7) 


(56) 


a) = 


1^7 7*^(7) 


(57) 


a^a = 




(58) 



it follows that the displacement operators form a complete set of operators 
for that mode. It is easy to generalize this proof to the multi-mode case. The 
displacement operators are over-complete. 



5 Operators 

Some operators can be written as sums of products of even numbers of cre- 
ation and annihilation operators; we shall call such operators even. Operators 
that can be written as sums of products of odd numbers of creation and an- 
nihilation operatorswe shall shall speak of as odd. Although most operators 
are neither even nor odd, the operators of physical interest are either even or 
odd. The number operator a^a, for example, is even, while the creation and 
annihilation operators, a' and a, are odd. 

The operators of quantum mechanics and of quantum field theory do not 
themselves involve Grassmann variables. Thus even operators commute with 
Grassmann variables, while odd ones anti-commute. 
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5.1 The Identity Operator 



If we compare the integral 



J d 2 a \a)(a\P) = J d 2 ct exp [a\a, + a* fa + a^a* + §A/3*) j |0> 

(59) 

with the integral-formula (|47|) and identify (3* and 7 in that formula with 
and j3 in this integral, then we have 

J d 2 cx |a)(a|/3) = exp ^ (at/3, + j |0) = 1/3). (60) 

Since the coherent states form a complete set of states, as shown by the 
expansion (]54j), it follows that the identity operator is given by the integral 

d 2 a\a)(a\. (61) 



The corresponding expression for the identity operator in terms of the 
eigenstates |a)' of the creation operators is 

/ = jUi-d 2 ^) \a)'\a\. (62) 



5.2 The Trace 



The trace of an arbitrary operator B is the sum of the diagonal matrix 
elements of B in the n-quantum states, 

Trfi = ^2(n\B\n), (63) 

n 

which shows that the trace of an operator that is odd vanishes. By inserting 
the preceding formula fl6"T|) for the identity operator, we have 

Trfi = I d 2 a{n\a){a\B\n). (64) 

n J 
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If we move the coherent-state matrix element (n\a) to the right of the matrix 
element of the even operator B, then we see from the formula fl22|) that 
minus signs arise that can be absorbed into the argument of either of the two 
coherent states, 

TrB = / d 2 a (a\B\n)(n\ — a) = / d 2 a (—a\B\n)(n\a), (65) 

n n 

in which the sum ^ n |^)(^| = / is the identity operator. The resulting 
multi-mode trace formula is 

TrB = Jd 2 a(a\B\- a) = j d 2 a (-a\B\a), (66) 

which holds also for odd operators, both sides vanishing. An important 
example is the trace of the dyadic operator |/3)( 7 |, 

Tr|/3)( 7 | = jd 2 cx{cx\l3){ 1 \-cx) = jd 2 cx{ 1 \-cx){ ( x\f3) 

= |rf 2 a(- 7 |a)(a|/3) = (- 7 |/3) = (7|-/3), (67) 

in which we have used the completeness relation Q5"T|). Since the coherent 
states are complete, we may replace in this formula either the ket |/3) or 
the bra (7I with its image F\/3) or (~f\F under the action of the arbitrary 
operator F and obtain the trace formula 

Tr(F|/3)( 7 |) = Tr(|/3)( 7 |F) = (- 7 |F|/3> = ( 7 |F| -(3). (68) 



5.3 Physical States and Operators 

A state is physical if it changes at most by a phase when subjected to a 
rotation of angle 2n about any axis, 

U(h,27r)\t/j) = e ie \^j). (69) 

Since fermions carry half-odd-integer spin, a state of one fermion or of any 
odd number of fermions changes by the phase factor —1. States that contain 
no fermions or only even numbers of fermions are invariant under such 2ir 
rotations. 
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Thus physical states are linear combinations of states with odd numbers 
of fermions or linear combinations of states with even numbers of fermions. 
But a state that is a linear combination of a state that contains an odd 
number of fermions and another that contains an even number of fermions 
is excluded. For instance, the state 

^(|0> + |1» (70) 

is unphysical because under a 2n rotation it changes into a different state: 

U(n,2n)^= (|0> + |1» = ±= (|0> - |1» + e*-L (|0> + |1» . (71) 

We define an operator as physical if it maps physical states onto physical 
states. Physical operators are either even or odd. 

In all physical contexts that have been explored experimentally, the num- 
ber of fermions (or more generally the number of fermions minus the num- 
ber of anti-fermions) is strictly conserved. That conservation law leads to 
certain further restrictions on the permissible states of the field. If we let 
N = J2 k a\dk be the fermion number, the law requires that any state arising 
from an eigenstate of N must remain an eigenstate of N. This law can be 
derived from an assumed U(l) invariance of all the interactions under the 
transformation U{9) = exp (iON), which changes a and a) to 

e -i9N ae iON = e i8 a (?2) 

and 

e -»N a 1 e iON = e -i6 a \ (73) 

Fermion conserving interactions involving the a& and a k are ones in which the 
phase factors e ±t6 all cancel. If a system begins in a state with a fixed number 
of fermions, the conservation law restricts the set of accessible states consid- 
erably more than the In super-selection rule mentioned earlier. Transitions 
can not be made, for example, between states with different even fermion 
numbers or between states with different odd fermion numbers. 
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5.4 Physical Density Operators 



A physical density operator can be written as a sum of dyadics of physical 
states with positive coefficients that add up to unity. It follows that a physical 
density operator p is a positive hermitian operator of unit trace: for any state 



W) > (74) 
P f = P (75) 
Tip = 1. (76) 



Physical density operators are invariant under a 2n rotation. Thus the 
one-mode operator 

P = ^(|0)(0| + |1)(1|), (77) 
for example, is a physical density operator, but the dyadic 

^(|0) + |1))((0| + (1|) (78) 

is not. In this work we shall consider only density operators that are physical 
in this sense. 

The dynamical problems we solve do not always begin with a fixed number 
of fermions. More generally they begin with a mixture of states with different 
fermion numbers, that is with density operators of the form 

P = J2pn>\N > )(N>\ (79) 

N' 

where the pw are real and non-negative. Such density operators are invariant 
under the transformation U (9) = exp(i9N), and the fermion conservation law 
assures us that they will always remain so, 

e- im pe im = p. (80) 



The coherent states do undergo a simple change under this transforma- 
tion, 

U(9)\a) = e ieN \a) = \e t6 a) (81) 
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(a\W(6) = (a\e- ieN = (e ie a\, (82) 
which leaves their scalar product invariant, 

(e ie a\e ie a) = (a\a). (83) 

6 Delta Functions and Fourier Transforms 

We can define a function 

8(t-C) = Jd 2 ccex V ^(a n (C-C)-(L-Cn)<)\ (84) 

= U^-Cn) (C-C) (85) 

n 

which plays the role of a Dirac delta function in that if /(£) is any function 
of the set £ of Grassmann variables {£i, £ 2) • • •}, then 

0/(0= /(C)- (86) 

The delta function is doubly even: it commutes with Grassmann numbers 
and*(e-C) = *(C-fl- 

We have been using the term Fourier transform to denote an integral of 
the form 

/(«) = Jd^e^-^fiO- (87) 



The delta-function identity (|84j ) implies that the inverse Fourier transform is 
given by the similar formula 

/(0 = Jd 2 ae^- a ^f(a). (88) 

The identity flS4| ) also leads to two forms of Parseval's relation: 

J d 2 a /(a) [g{a)T = J /(0 <?*(0 (89) 
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and 



J d 2 a /(a) g(-a) = J d 2 £ /(£) <?(£), (90) 

which apply also to operator-valued functions provided that complex conju- 
gation is replaced by hermitian conjugation. 



We may use the formula (84) for the delta function to derive a fermionic 



analog of the convolution theorem: 

d 2 f3f(a-(3)~g((3), (91) 



which expresses the Fourier transform of the product of the two functions 
/(£) and g(£) as the convolution of their Fourier transforms f(a — f3) and 

By using the normally ordered form ([16]) of the displacement operator, 
the eigenvalue property of the coherent states, and the preceding formula 
(|84T) for the delta function we find 



d 2 7 (7|£>(a)|7> = /c/ 2 7(7|e ata e- a * a | 7 )e2 c 

d 2 J e 7*a-a*7+^«"* 



= 5(a)e* aa =5(a). (92) 

The addition rule ([19]) for successive displacements now implies that for the 
multi-mode case 

J d 2 ~f ( 7 \D( a )D(-f3)h) = 8(a - (3). (93) 



7 Operator Expansions 



The preceding delta-function identity fl93|) and the completeness (|55]-|58|) of 
the displacement operators give us a means of expanding an arbitrary oper- 
ator F in the form 

F = Jd 2 if{i)D{-i). (94) 
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We may solve for the weight function /(£) by multiplying on the right by 
the displacement operator D(ot) and then taking the diagonal coherent-state 
matrix element in the state \/3) and integrating over (3: 



d 2 (3(f3\FD( a )\f3) = jd 2 (5 Jd 2 tm(P\D(-t)D( a )\(3) 

= Jd 2 Zf(£)5( a -£) = f(cx). (95) 

The full expansion is thus 

F = jd 2 i Jd 2 f3((3\FD(t)\f3}D(-Z). (96) 
Such expansions will prove useful in the sections that follow. 



The formula (p4|) for the delta function <5(£ — () may be interpreted as a 
trace identity. From the eigenvalue property of the coherent states, it follows 
that 

*<{ - C) = = /^e— W^V^a) (97) 

in which we recognize the normally ordered form ( |T5] ) of the displacement 
operator 

S(£ ~ C) = J d 2 ae a ^ a * (a\D N (()\a). (98) 
By using the trace formula (0), we may write this delta function as the trace 

6(Z-C) = Jd 2 ae a ^-^Tr[D N (()\a)(-a\] 

= Tr[D N {()E A (-Z)] (99) 

The subscript A has been chosen to indicate its anticipated of the product 
of the normally ordered displacement operator Dn(Q with an even operator 
Ea{0 defined as the Fourier transform 



(100) 



of the coherent-state dyadic |a}(— at\. As intimated by its subscript, the 
operator E A (^) will turn out to be useful for dealing with anti-normally 
ordered operators. 
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We may now use the completeness (pHHoTf) of the displacement operators 
and the trace identity (|99|) to expand an arbitrary operator F in terms of 
the normally ordered displacement operators Dn(£), 

F = Jd 2 U(0D N (-0- (101) 

We may solve for the function /(£) by multiplying on the right by the oper- 
ator E A (() an d forming the trace: 

Tr [FE A (0] = J d\ /(OTr [D N {-^)E A {Q] = J d\ f{t)5{C - = /(C)- 

(102) 

The full expansion is thus 

F= ( d 2 |Tr [FE A (t)] D N {-£). (103) 



By using the Grassmann calculus, one may compute the Fourier transform 
( lOOp of the coherent-state dyadic |ck)(— oc\ and find for the operator E A (£) 
the formulas 

E A (0 = |0)(0|-(r + a t )|0)(0|(^ + a) (104) 
= 2(i - da) + £fW + fa* - £*a (105) 

with which it is easy to exhibit the completeness of the operators E A (^): 

d 2 i 2(1 + C0E A (0 (106) 
d^(-0E A (0 (107) 
at = jd^(-C)E A (0 (108) 
i-ota = jd 2 i\iCE A {i). (109) 



Since the operators E A (£) are complete, we may expand an arbitrary 
operator G in terms of them, 

G = jd^g(0E A (-0 (110) 
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and then use the trace formula (P§| ) and the evenness of the displacement 
operators to evaluate the weight function 

Tr[D N (()G} = J<Ptg{t)Tr[D N {QE A (-£)] = J <?tg(t) 8(£ - Q = g(C) 

The full expansion is thus 



1111 



G = Jd 2 $Tr[GD N {£)) E A {-£). 



(112) 



8 Characteristic Functions 



For a system described by the density operator p, we define the characteristic 
function x(£) °f Grassmann argument £ ( and £* ) as the mean value 



X(t) = Tr 



;ii3) 



It is thus a species of Fourier transform of the density operator p. Because 
£? = £i* 2 = 0, we may expand the exponential as 



X(0 = Tr 



P II ( X + ~ a ™^« + £&i(4°n - I)) 



;ii4) 



We may also define the normally ordered characteristic function Xn{£,) 



as 



Xn{£) = Tr 
with the expansion 

Xn{£) = Tr 



p exp £„4 exp 



PYI( 1+ ^ nCL n ~ a nCn + C£n4°n) 



;ii5) 



(116) 
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The anti-normally ordered characteristic function Xa{£) is 



Tr 



p exp 



exp 



117) 



Tr 



I"J (1 + £n4 - a n £* + C^(4a„ - 1)) 



118) 



Because the density operator p is an even operator and because the displace- 
ment operators are constructed from bilinear forms in fermionic quantities, 
it follows that the characteristic functions are even in the sense that they 
commute with Grassmann variables. 



8.1 The S-Ordered Characteristic Function 

We may define a more general ordering of the annihilation operator a n and 
the creation operator a*, much as we did earlier for boson-field operators 
It is an ordering specified by a real parameter s that runs from s — — 1 for 
anti-normal ordering to s = 1 for normal ordering. For the quadratic case, 
the s-ordered product for fermions is 

{ala n } s = ala n + ~(s - 1), (119) 

to which we append the trivial definitions 

{al} s = a\ and {a n } s = a n . (120) 

We note that the definition ( |119| ) differs by a crucial sign from that 0] of 
s-ordering for bosonic operators b n and b^: 

{bib n } s = bib n + \{i~s). (121) 

In particular the anti-normally ordered product {a^a n } and 
the symmetrically ordered product {a^a n } is half the commutator, 

W n a n } = \ [a n , a^] . (122) 
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We define the s-ordered characteristic function s) as 



Tr 



Tr 



Tr 



p II i 1 + ~ a ™^« + W n a n } s ) 

p n ( x + - + ( a « a « + a* - !))) 

P exp I ^ ((„4 - a«C + | C£n 

\ n 



(123) 
(124) 
(125) 
(126) 
(127) 



which, incidentally, shows it to be an even function. 



A particularly useful example of these characteristic functions is the case 
of the anti-normally ordered function xa(£) = x(£> ~ !)• We see by inserting 
the resolution (|6lD of the identity between the exponential functions in its 



definition (|117 ) that 



*(£,-!) = Tr 



p exp I - PmC m J J d 2 (3 \(3)(f3\ exp ( ^ £ n /£ 



128) 

in which we have replaced the annihilation and creation operators by their 
eigenvalues in the coherent states. By using the trace formula (|57D, we find 



(129) 



which expresses the anti-normally ordered characteristic function x(£, —1) as 
the Fourier transform of the matrix element ((3\p\ — (3). 



If we define the s-ordered displacement operator D(£, s) as 



(130) 
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then we may write the s-ordered characteristic function (|127|) as the trace 

X (t,s) = Tr[pD{t,s)]. (131) 



9 S-Ordered Expansions for Operators 

A convenient extension of the definition of the operator E A (£) is 

E{£, s) = E A (£) exp (^±± £ ejr) , (132) 
from which we note that 

E A (t)=E(£,-l). (133) 

This is one sense in which the operator E A (£) is related to anti-normal or- 
dering. 

By using the s-ordered operators D(£,s) and E(£,s), we may generalize 
the expansions ( p. 031 ) and (|112|) of the arbitrary operators F and G to 



F = Jd 2 ZTr[FE(Z,-s)]D(-£,s) (134) 

G = Jd 2 £Tr[GD{€,-s)} E(-£,s). (135) 

The obvious generalization 

8{t-C) = Tr[D(£,8)E{-C-8)} (136) 
of the trace formula (SDp then gives the trace of the product FG as 

Tr [F G] = Jd 2 dTr[FE(t,-s)] Tr [GD(—£, s)] . (137) 

We may now use the second Parseval relation (|90|) to cast the expansions 
(1134j ) and ( |135| ) into forms that will prove to be quite useful. First let us 
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define the complete sets of operators D(cc,s) and E(a,s) as the Fourier 
transforms of the operators D(£, s) and E(£, s): 

D(a, s) = J d 2 £ exp (a n £ - £n<) j s) (138) 
= j d 2 d exp ^ (a n £ - Cn<) j s). (139) 



Ma, s 



Next let us define the weight functions F E (oc, —s) and Gd(oc, —s) as the 
Fourier transforms of the traces 



/ d2 z ex p (e - j Tr t Fi ^> 
/ rf2 £ ex p (E - Tr • ( 141 ) 



= / <l 2 £ exp ( X (m„C - CO ) Tr [/•'£'(£ -«)] (110) 

Gd(ql, -s) 



It follows then from the Parseval relation flDTf ) and from the expansions 
and (|135|) that the operators D(ct,s) and E(oc,s) form complete sets of 
operators and afford us the expansions 



draF E (a,-s)D(a,s) (142) 
G = I d 2 ctG D (a, -s) E(a, s) (143) 



of the arbitrary operators F and G. Applying the Parseval relation fl9"0| ) to 
the trace formula ( 137|) , we have the trace relation 



Tr[FG]= d 2 aF E (cx,-s)G D (ct,s). (144) 



The operators E(a,s) are particularly simple when s = ±1. It follows 
from the definitions (|132|) and Q100Q of the operators E(oc, s) and Ea(£), and 
from the formula fl8"4] ) for the delta function that the operator E(oc, —1) is 
just the coherent-state dyadic 

E(a,-1) = \a)(-a\. (145) 
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Similarly, by using the definitions ( |132| ) and ( |10(J| ) and the Fourier-transform 
relation fl48|) , one may write the operator E(a, 1) as the integral 

E( a , -1) = /J! H 2 &) e-l a -P>< a *-F) \(3}(-f3\. (146) 

By performing the integration and referring to the explicit formula ([H]), we 
may show that the operator E(a, 1) is the dyadic of the eigenstates (p7| ) of 
the creation operators |a)': 

£(a,l) = |a) M (-a|. (147) 



10 Quasi-Probability Distributions 

Among the most important of the foregoing expansions, is the expansion 
(|135|) when the operator G is the density operator p, 



p = Jd 2 £Tr[pD(Z,s)]E(-Z,-s), (148) 
in which case the trace is the s-ordered characteristic function x(£i s )> 

P= fd 2 Zx(ts)E(-Z,-s). (149) 



We may define the s-ordered quasi-probability distribution W(oc, s) as 
the Fourier transform of the s-ordered characteristic function x{& s ) 

W{a,s) = ^exp ^KC-CnOJ X(t,s); (150) 

both W(ot, s) and x(^, s) are even functions. It follows now from the expan- 
sion ( p_43| ) that the s-ordered quasi-probability distribution W(a, s) is the 
weight function for the density operator p in the expansion 

p = Jd 2 otW(a,s)E(ot,-s). (151) 
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The functions W(a, s) for different values of the order parameter s are 
intimately related to one another because the characteristic functions obey 
the identity 

(a -=^r-€) (152) 



X (t s) = exp ( -£* • £ ) x(£) = exp 



The function W(ct, s) is therefore the Fourier transform of the product of 
exp ^^^^* • with the characteristic function x(£>t) 



W(oc, s) = J d 2 $, exp KC - Cn<) j exp ^ 



(153) 

The Fourier transform of the characteristic function x(£>^) is ^(c*^); while 
that of exp • ^) according to Eq.® is 

-(t-s) 



n 

The convolution theorem (PTJ) now gives s) as 

W(a,a) = |n 



■ e(*- s ) 



7"7n 



(t- S ) 



(t- S ) 

A useful example of W(a, s) is the function 



(154) 

W((3,t). 
(155) 



W(a,-1) = ^exp ^KC-a<)j X(^-1), (156) 



which according to Eq.(129) is the Fourier transform 



W(a,-1) = J d 2 £d 2 /3 exp 



((3\p\-(3). 
(157) 

By using the delta-function identity (|84j), we see that this expression reduces 
to 



W{a,-1) = Jd 2 f35( ( x-f3)((3\p\-(3} 
= (ct\p\ — a). 



(158) 
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This function is the fermionic analogue of the function Q(f3) = (f3\p\f3) which 
is often used to represent the density operator p in terms of the bosonic 
coherent states |/3). It is the weight function that gives the mean values 
of anti-normally ordered products of creation and annihilation operators in 
terms of integrals of the corresponding products of Grassmann numbers. 



11 Mean Values of Operators 



We shall here be concerned with computing the mean values of the products 
of s-ordered monomials 

/ . \ n ; 

(159) 



in which the exponents rii and take the values or 1. The ordering of 
the modes labelled by the index i is arbitrary but fixed. We shall show that 
we may express the mean values of such products of monomials as integrals 
of the s-ordered weight function W(ct, s) multiplied by the monomials in the 
same order. By using the definition ( |150| ) of W(ct,s), we may write these 
integrals in the form 

/ d 2 cx J] (a*r «r W{a, s) 

= Jd 2 cc\{ (a*) ni aT J d 2 £ex P ^ - J x (t, *).(160) 

It is now easy to write the monomial as a multiple derivative, 
J d 2 cx J] (a*) ni a? W(a, s) 



-e * ! 1 



On using our formula ([52]) for integration by parts, we have 



X(ts). (161) 



a 



d 2 £d 2 aexp ^ (a^* - ^a*) ] [ 



*\m,i 
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in which we recognize the delta-function formula fl55]) which gives 



jd 2 cx\[{a*rarW{ a ,s) 



2 Q n i Q m i 

d ££(£) H Q(jTyH d(Q) mi 



11 Q^yT, Q^rrn X ^ ^ 



TV 



(163) 
(164) 

(165) 



If we recall the definitions of s-ordering in Eqs.( |119|) and ( 120 ), we then find 



J d 2 a J] (a*r «r W(a, s) = Tr p \{ { ^ a?*} 



(166) 



In particular, by taking n,i = rrii = 0, we see that the weight function 
W(ot, s) is normalized, 



W(ct,s) = Tip = 1. 



(167) 



12 The P-Representation 



Of the representations ( |151| ) for the density operator p, by far the most 
important is the one for s = 1 with the normally ordered weight function 
P(ct) = W(ct, 1). By Eq. ( |145| ) it takes the simple form 



p= d 2 at P(ct) \ct)(—ct\, 



(168) 



which recalls the P representation for boson fields. Since the function 

P(ct) is even, we may also write 



p— j d 2 ct P{ol) I — a) (oc\. 



(169) 
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Because Grassmann integration is differentiation, the fermionic P represen- 
tation is not affected by the mathematical limitations ||— 1| that restricted 
somewhat the use of the bosonic P representation. 

The P-representation may be used directly to compute the mean values 
of normally ordered products 

Tr^pafraf) = j d 2 cxP(cx){cx\a^a[ n \a) 



d 2 aP(a)a>af- 



(170) 



This extremely useful relation is just a special case of Eq.( |166| ) for s = 1. 



Since the operator E(ac, 1) is the dyadic ( |147| ) of the eigenstates of the cre- 
ation operators, it follows from the expansion ( |151| ) that the weight function 

flU) ' 

Q{a)=W{a,-l) = (a\p\-a) (171) 



is the weight function in the representation 



P 



d 2 ctQ(cx) |a) M (-«l 



(172) 



which affords the simple way of computing the mean values of anti-normally 
ordered products that corresponds to Eq. (|166|) for s — — 1. 

Another use of the weight function Q(a) = W(ct, — 1) = (oc\p\ — at), 
however, is that it allows us to compute the weight function P(cx) = W(a, 1) 
of the P-representation as the simple convolution 



P(a) = f\{ {-d 2 (3 m ) exp 



((3\p\-(3), (173) 



as follows from the general convolution formula (|155|) with s — 1 and t — — 1. 
Although the analogous relation for bosons often is singular 0-0, this result 
holds for all fermionic density operators p. 
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13 Correlation Functions for Fermions 



A principal use of the P representation for bosonic fields has been the evalu- 
ation of the normally ordered correlation functions, which play an important 
role in the theory of coherence and of the statistics of photon-counting experi- 
ments 0. The analogously defined correlation functions for fields of fermionic 
atoms can be shown to play a similar role in the description of atom-counting 
experiments 0. If we use i/j(x) to denote the positive-frequency part of the 
Fermi field as a function of a space-time variable x, then the first two of these 
correlation functions may be defined as 

G^(x,y) = Tr [^(x)^(y)] (174) 

G {2 \x 1 ,X2,y 2 ,yi) = Tr [pi/j ] (x 1 )^(x 2 )^(y 2 )i>(yi)] • (175) 

(176) 

The nth-order correlation function is 
G (n) (xi, ...,x n ,y n ,...,y{) = Tr [p^fai) • • • ^(x n )if)(y n ) . . . 4>(yi)] ■ (177) 

If we expand the positive-frequency part of the Fermi field in terms of its 
mode functions (f>k(x) as 

4j(x) = J2 a kMx), (178) 

k 

then its eigenvalue in the coherent state \ac) 

i(){x)\a) = (p(x)\a) (179) 

is the Grassmann field 

<p(x) =J2 a Mx) (180) 
k 

in which the annihilation operators in ( |1 78| ) are replaced by the Grassmann 
variables a = {a k }. 

We may use the P representation to evaluate the nth-order correlation 
function as the integral 

G<> n \x 1 ,...,x n ,y n ,...,y 1 ) (181) 
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= Jd 2 ccP{a){oL\^\x l )...^\x n )ij{y n )...^{y l )\oL) (182) 
= Jd 2 OLP(a) if ] {xi) . . . ip ] (x n )(p{y n ) . . .(p{yi). (183) 

14 Chaotic States of the Fermion Field 

The reduced density operator for a single mode of the fermion field can be 
represented by a 2 x 2 matrix for the states with occupation numbers and 
1. If the matrix is diagonal, it is specified completely by the mean number 
of quanta (n) in the mode. The density operator for the kth mode, in other 
words, must take the form 

Pk = (l-(n k ))\0)(0\ + (n k )\l)(l\. (184) 

We shall speak of this density operator as representing a chaotic state of the 
kth mode. A chaotic state of the entire field will then be represented as a 
direct product of such density operators for all the modes of the field, 

Pch = JJpfc- (185) 

k 

It is specified by the complete set of mean occupation numbers {(^fc)}. 

The total number of fermions, N = a\ak, present in chaotic states will 
in general be indefinite. Indeed it is easily seen that in the state specified by 
Eq. ( |185| ) we have 

<^ 2 >-<A0 2 = 5> fc )(l-K>) (186) 

k 

so that N can not be fixed unless all the (n^) take the values or 1. The 
indefiniteness of the number of particles present is a feature that the chaotic 
states of the fermion and boson fields have in common. For sufficiently large 
values of N, however, the fluctuations of N/ (N) may be quite small so the 
specification of N in these relative terms may be quite precise. Fluctuations 
of this type in the number of particles present are a familiar property of the 
grand canonical ensemble in statistical mechanics, and that ensemble, as we 
shall see, represents a special class of chaotic states. 
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The single-mode density operator ( |184| ) can also be written as 

Pk = (1 - K» (t^TtX^ C|0><0| + |1)(1|) (187) 



1 - (n k ) 

in which we recognize the unit operator I k for the subspace of the kth mode. 
Within this subspace we have 

ft = (1 - W) (i^))" ! *' (188) 

This expression can be used quite directly to evaluate the weight function 
W{a,-l) = Q{a). 

We first note that for any real number v 

v ata \a) = e^ aa * {1 ~ v2) \av), (189) 



so that we have 



(a\v aia \ -a) = e V*(i-* 2 ) ( a \ - av) 

= e aa ' {1+v) . (190) 



Then if we let v = (rife) /(l — (n k )), we see that 

1 - <"A-) 

and 



(a k \Pk\ ~ a k ) = (1 - (n fc )) exp ( - ^ - ) (191) 



Q(a) = W(a,-l) = Y[(a k \p k \- a k ) 

k 

- lid - exp (^_) . (192) 

This product is the weight function appropriate to averaging anti-normally 
ordered operator products in chaotic states. 

We may find the weight functions corresponding to all the other ordering 
schemes by using the convolution ( 155 ) with t = —1 and carrying out the 
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required integration with sufficient attention to the implicit minus signs. The 
result for the kth mode is 

s + 2{n k ) - 1 ____ ( 2a k a* k 



Wk(at , .) = exp ^ s + 2( ;;_ 1 ) , (193) 

and the weight function for the multi-mode field is simply the product 

W(cx,s) = l[W k (a k ,s). (194) 



Thus the function W k (a k ,0), which is analogous to the Wigner function 
for boson fields, is given by 

W k (a k , 0) = - «n fc > - |) exp {-^^j , (195) 

and the function W k (a k , 1), which is the analogue of the function P k (a k ) for 
boson fields, is 

W k (a k , 1) = P k (a k ) = -(n k ) exp ("^y) ■ (196) 

The latter result is a particularly useful one since there are many physical con- 
texts that call for the averaging of normally ordered products of annihilation 
and creation operators. For chaotic fields one may calculate all such averages 
as Grassmann integrals by making use of the fermionic P-representation with 
P(a) given by Eq.( |19g ). 



The minus signs in front of the expressions ( |195| ) and ( |196| ) may be some 



what surprising since these functions are the fermionic analogues of quasi- 
probability densities that are predominantly positive for boson fields. It is 
worth pointing out, therefore, that these signs result from our convention 
that defines d 2 a as da*da. Had we chosen the differential instead to be 
dada*, the signs would have been positive. 

For a chaotically excited boson field, the P-representation expresses the 
density operator as a gaussian integral of a diagonal coherent-state dyadic. 
For fermion fields the corresponding expression of p k for a single mode is 

f a k a % 

Pk = -(n k ) d 2 a k e~^ \a k )(-a k \. (197) 
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According to Eq.( 187 ), the density operator p k can also be written as a sum 
over the m-fermion states as 

Pk = (l-(n k )) J2 ) ) k K>Ki- (198) 

What we have shown, in effect, is that the two expressions are identical and 
that statistical averages can be evaluated by means of gaussian integrations 
for fermions as well as for bosons. The multi-mode density operator is rep- 
resented, of course, by the product of the single-mode density operators, 

P = UkPk- 

Fields in thermal equilibrium with a suitable particle reservoir represent 
particular examples of the kind of chaotic excitation we have been describing. 
If it is appropriate to describe such fields by means of the grand canonical 
ensemble, then their overall density operator may be written as 

P=^—,e-^ N \ (199) 

where f3 = l/ksT, p is the chemical potential, H is the hamiltonian for the 
system, N is the particle number, and the normalizing factor H(/3, p) is the 
grand partition function. For a field with dynamically independent mode 
functions labeled by the index k, we can write 

H = ^e k a\a k , N=J2 a W (200) 

k k 

where e k is the energy of a particle in the fcth mode. 

Under these circumstances the equilibrium number of fermions in the fcth 
mode is 

<"*> = < 201 > 

In that case the ratio (n k )/(l — (n k )) is simply the generalized Boltzmann 
factor 



= e 0to-*O. (202) 
1 - (rife) 

We then find that the product of the p k given by Eq. (|197|) is precisely equal 
to the grand canonical density operator 



J] (-(n k )d 2 a k e-^T^j \a)(-a\ = e'^ H ^ N \ (203) 
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There are many examples of thermal equilibria for which the P-representation 
on the left should furnish a useful computational tool. 



15 Correlation Functions for Chaotic Field 
Excitations 



We have introduced a succession of normally ordered correlation functions 

G (n) 

(xi . . . x n , y n . . . yx) in section 13, and shown how they can be expressed 
as integrals over the Grassmann variables ct = {a k }. For the case of chaotic 
fields, the appropriate weight function is 

P{a) =Y[P k {a k ), (204) 

k 

the product of the gaussian functions in Eq.( |196j ). The first-order correlation 
function is thus given by 

G (1 \x,y) = Jl[^-(n k )d 2 a k e-^y^ (a\^{x)^{y)\a). (205) 

The fields ip and ijy may now be replaced by their Grassmann field eigenvalues 
defined by (|179|) and ( |180|) . Their product is a quadratic form in the variables 



a k and al, which is easily integrated: 

JJ f -(n k )d 2 a k e'~ J ^ a i^m<f>j(x)(f> m (y) 

k l,m 

= ^2(n k )<l>l(x)My)- ( 206 ) 

k 

To find the higher-order correlation functions, we can make use of a 
species of generating functional. We first define the Grassmann fields 

C(x) = ^/^(x) (207) 
k 

rj(y) = $>^fc(y), (208) 
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and use them to construct the normally ordered expectation value 



r[c,77] = iv 



pexp ^ C(x)ip i (x)d 4 x S j exp (^j ij(y)r]*(y)d 4 y 



(209) 



If we form the variational derivative of T with respect to C(^i) from the left 
and with respect to rj*(yi) from the right, subsequently setting ( and i] to 
zero, then we find an alternative expression for the first-order correlation 
function, 

sik) iJm r| <— = Tr ^'w^a = G<1,( *"^ (210) 

where left and right differentiation have been indicated explicitly in the sub- 
scripts. 

It is evident then that one may generate all of the higher-order correlation 
functions by performing further differentiations, 

r ( n ) ( s _ 5 5 5 5 

[Xl /7 fcC(si) ' ' ' hCM S RV *(y n ) ■ ■ ■ 5 RV *( yi ) '^=°- 

(211) 

To evaluate the generating functional T for a chaotic field, we make use of the 
orthonormality of the mode functions (pk and then carry out the Grassmann 
integration 

r = j\{ ^-(n k )d 2 a k e~^y^j exp + atf) 

= II + K)^7fc) = exp I ^{njPkrt 

= exp C^G^y^d 4 ^) • (212) 

If we begin performing the variational differentiations to find the second-order 
correlation function, we may write 

5 5 



C(x)G^(x,y 2 )d 4 x I ((x')GW(x', yi )d 4 x'. 

(213) 
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We then find 



G (2) (x 1 x 2 y 2 yi) 



5 5 




and since ((x) and C(x') anti-commute, 



G {2) (x l ,x 2 , y 2 , Vl ) = G^(xi, yi )G<U(x 2 , y 2 ) - G^(x u y 2 )G^(x 2 , Vl ). (215) 



The generalization to nth order is immediate. It expresses the nth-order 
correlation function for chaotic fields as a sum of products of first-order 
correlation functions with permuted arguments, 



This expression is summed over the n! permutations of the indices 1 . . . n. 
The factor (— l) p is the parity of the permutation, and the index Pj is the 
index that replaces j in the permutation. 

The expression of the nth-order correlation function in terms of first- 
order correlation functions is characteristic of chaotic fields. Such fields are 
completely specified by the set of mean occupation numbers (n k ), and these 
are already contained in the first-order correlation function. 



16 Fermion Counting Experiments 

The use of photon counting techniques has for many years been the most 
direct means of investigating the statistical properties of light beams. Ex- 
periments of this type began with that of Hanbury Brown and Twiss || 
in 1956 and expanded greatly in scope with the development of the laser. 
The theory || underlying these experiments is based on the evaluation of 
quantum-mechanical expectation values of normally ordered products of elec- 
tromagnetic field operators. The coherent states of the field thus play a 
special role in the formulation of that theory. The application of the theory, 
furthermore, extends to boson fields of much more general sorts, including 
for example beams of heavy atoms |J. 



n 



G^\x x ...x n ,y n ... Vl ) = £(-l) p J] G{1 \*v VP,)- 



(216) 



P 3=1 
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In the case of the electromagnetic field, it has been shown || that the 
probability of detecting n photons in a given interval of time can be expressed 
as the nth derivative with respect to a parameter A of a certain generating 
function Q(A), 

(—l) n d n 

M-hrw™"- (217) 

The generating function Q(X) for the electromagnetic field is the expectation 
value of a normally ordered exponential function of the form 

Q(A) = Tr (p:e- AJ :) , (218) 

in which the symbols : : stand for normal ordering, and the operator I is 
a space-time integral of the product of the positive-frequency and negative- 
frequency parts of the field, and E^~\ respectively. 

For the case of fermion fields, it can easily be shown || that the prob- 
ability of counting n fermions in a given interval of time falls into precisely 
the same general form. In the simplest instance, for detectors that respond 
to the density rather than the flux of the particles, the integral 1 takes the 
form 

J = K ^(r, t)ij){T, t)d 3 rdt } (219) 



where the constant k is a measure of the sensitivity of the counter and the 
integration is carried out over the counting-time interval and over the volume 
being observed. 

To obtain the expectation value of the normally ordered exponential func- 
tion in Eq.( |218| ), then we may use the P- representation for the density opera- 
tor p. In that case the field operators ijj(f, t) and ^(f, t) are, in effect, always 
applied to their eigenstates, coherent states such as \cx) and (a\. They can 
then be replaced by their Grassmann field eigenvalue functions defined by 



Eq.( |180| ) and its adjoint, so that we have 

Q(A) = / d 2 ccP(a)e- XJ , (220) 



where 

J = k I cp*{r, t)cp(r, t)d 6 rdt. (221) 
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The expression J is a quadratic form that we can write as 

J = J2 a *k B kk' a k'> ( 222 ) 



k.k' 



so the evaluation of the generating function Q(A) reduces to the calculation 
of the integral 

Q(A) = ( d 2 aP(a) exp ( -A ^ a* k B w a k , ) , (223) 

in which the normal ordering symbols are no longer necessary because of the 
simple anti-commutation properties of the Grassmann variables a k . 

For the case of the chaotic fields defined in section 14, this integral takes 
the form 

Q(A) = /il {-(n k )d 2 a k e^\ exp I -A £ a* k B ky a v ) . (224) 

If we define a new set of variables (3 k = a k j \J (n k ), we find according to the 
rule ( f4"o] ) that the integral can be written as 

Q(A) = / l[(-d 2 (3 k ) exp ((J fcfc , - \M kk ,){3 k \ , (225) 

fe \k,k> J 

where the matrix M is 

M kk i = (n k )B kk ^ (n k >). (226) 

A unitary linear transformation on the variables f3 k can then be used to 
diagonalize the quadratic form in brackets. If the eigenvalues of the matrix 
1 — AM are /i;, then the integral is easily seen, according to the formula (fHf ) 
for a = 0, to be 

Q(A) = Y[ in = det (1 - AM) . (227) 

This result may be used directly to find the various probabilities given by 
Eq. fl2171 ). It contrasts quite interestingly with the generating function for 
boson counting distributions, which with closely corresponding definitions 
takes the form 

= det (1+ AM) " (228) 
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17 Some Examples 



17.1 The Vacuum State 



For the density operator 



p= |0...0)(0...0| 



(229) 



which represents the multi-mode vacuum state, the normally ordered char- 
acteristic function Xn(£) is 



N 



(0 . . . 0| exp Cn4 j exp ^- £ a n & j |0 . . . 0) = 1. (230) 



The weight function of the P representation is then 



'(a) = J d 2 d exp ^ fas - &<) j = 5(a)- 



(231) 



The mean values of the normally ordered products of creation and annihila- 
tion operators all vanish 



Tr 



n aT ' =j d2a ii ^ r 6{a) = ° (232) 



except for the trace 



Tr[p] = J d 2 a5(a) = I. 
The general weight function W(a,s) of the vacuum is given by 



W(a,s) = -(1 -s) exp f 



(233) 



(234) 
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17.2 A Physical Two-Mode Density Operator 



Let us consider the most general physical two-mode fermionic density oper- 
ator 

p = r |00)(00| +u|10)(10| +u|01)(01| +w|10)(01| + w*|01)(10| 

+ x|00)(ll| + x* |11)(00| +t 1 1 1> < 1 1 1 (235) 

in which 1 10) = a{ |00), |11) = a^a^lOO), etc., and the Latin letters r,t,u, and 
v represent non-negative real numbers, while x and w may be complex. The 
non-zero traces are 

1 (236) 
(237) 
(238) 
(239) 
(240) 
(241) 
(242) 
(243) 

If the fermion number N commutes with the density operator p, then x = 
x* = 0. 

The normally ordered characteristic function Xn(£) is 

x(£)n = Tr P (i + £i a \ - a i£i + ti£i a \ a i) (l + 6<4 - °2^2 + Qtzoi *) 
= 1 + <*6 + ^^2^1 + (u + + (v + f)& 2 

+ x 66 + x* C2C1 + t Stefa. (244) 

The weight function W(oc, 1) is the Fourier transform of the normally 
ordered characteristic function Xn(£) 

W(a,l) = Jd^^il + a^ + a^i + ata^i) 

x (1 + a 2 £*2 + "26 + a* 2 a 2 ^) 

x[l+w £6 + w* &i + (u + t) Citi + (v + t) e 2 *6 



Tip 


= r + u + v + t 


Trpa\ai 


= u + t 


Tipa\a2 


= v + t 


Tipa\a\ 


= w 


Tr pa\ci2 


= w 


Trpaia2 


= X 


Tr pa\a\ 


= —x 


Trp a 2 a\ a\a 2 


= t. 
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and after following the rules (33-10), we find 

W(oc, 1) = t + wa 2 a\ + w*a\a 2 + (v + t)a\ai + [u + t)a 2 a 2 

+ x a±a 2 + x* a^al + a\a\a\a 2 . (245) 

We may now use this weight function to compute the mean values 
Jd 2 a 1 d 2 a 2 W(cx,l) 
Jd 2 aid 2 a2 a\a\ W(a, 1) : 
Jd 2 aid 2 a2 a 2 a 2 W(a, 1) - 
J d 2 aid 2 a 2 ol* 2 ol\ W(a, 1) 
Jd 2 aid 2 a 2 a\a 2 W(a, 1) : 

y ' S0L\d 2 0L 2 (X\OL 2 W [ol,X] ■ 

J d 2 a 1 d 2 a 2 ala 2 W(ct,l) ■■ 
d 2 a.\d 2 a 2 o^of x a.\a. 2 W(ot, 1) - 



Trp = 1 




(246) 


Trpa\ai 


= u + t 


(247) 


Tip a\a 2 


= v + t 


(248) 


Trpa 2 ai 


= w 


(249) 


Tr pa\a 2 


* 

= w 


(250) 


Tr pa\a 2 


* 

= X 


(251) 


Tr pa\a\ 


= —x 


(252) 


Tr pa\a\ 


a\a 2 = t 


(253) 



which agree with the results Q236| - |243"l) . 

With Pioi) = W(ct, 1) as given by ( |245| ) , we may write the density oper- 
ator ( |235| ) in the form of the fermionic P representation 



p= d 2 ccP(cx) \-a)(a\ = / d 2 et P(a) \a)(-a\. (254) 
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